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Abstract 

A study of Z-boson pair production in e + e~ annihilation at center-of-mass energies between 
190 GeV and 209 GeV is reported. Final states containing only leptons, and £ + £~vT'), 

quark and lepton pairs, (qq£ + l~ , qqVz/) and only hadrons (qqqq) are considered. In all states with at 
least one Z boson decaying hadronically, lifetime, lepton and event-shape tags are used to separate bb 
pairs from qq final states. Limits on anomalous ZZ7 and ZZZ couplings are derived from the measured 
cross sections and from event kinematics using an optimal observable method. Limits on low scale 
gravity with large extra dimensions are derived from the cross sections and their dependence on polar 
angle. 
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1 Introduction 



LEP operation at center-of-mass energies above the Z-pair production threshold has made a careful 
study of the process e + e~ — * ZZ possible. In the final data taking runs in 1999 and 2000 LEP delivered 
an integrated luminosity of more than 400 pb -1 to the OPAL experiment at energies between 190 GeV 
and 209 GeV producing of the order of 100 detected Z-pair events. In this paper we present results 
on Z-pair production from these data. Previously published studies from OPAL and the other LEP 
collaborations can be found in References [1-4]. 

In the Standard Model, the process e + e~ — ► ZZ occurs via the NC2 diagrams [5] shown in Figure 1. 
The Z-pair cross section depends on properties of the Z boson (the Z mass, raz, the Z resonance width, 
Tz, and the vector and axial vector coupling of the Z to electrons, <7y and g\) that have been measured 
with great precision at the Z resonance [6]. The expected Z-pair cross section increases from about 
0.25 pb at y/s = 183 GeV to slightly more than 1.0 pb above y/s = 200 GeV, but remains more than 
an order of magnitude smaller than the cross section for W-pair production. In contrast to W-pair 
production, where tree level WW7 and WWZ couplings are important, no tree-level ZZZ and ZZ7 
couplings are expected in the Standard Model. However, physics beyond the Standard Model could 
lead to effective couplings [7] which could then be observed as deviations from the Standard Model 
prediction in the measured Z-pair cross section and kinematic distributions. Such deviations have 
been proposed in the context of two-Higgs-doublet models [8] and in low scale gravity theories [9] . In 
this paper we report on measurements of the NC2 Z-pair cross section, including the extrapolation 
to final states where one or both Z bosons have invariant masses far from mi- These measurements, 
together with optimal observables determined from the final states of selected qqqq and qq£ + £~ events, 
are then used to extract limits on possible ZZZ and ZZ7 anomalous couplings. Finally we use the 
measured cross sections, binned in polar angle, to determine limits on large extra dimensions in low 
scale gravity theories. 

The analyses and presentation of this paper largely follow that of the first OPAL e + e~ — * ZZ 
results at center-of-mass energies of 183 GeV and 189 GeV [3]. All of the analyses have been optimized 
for the higher energies and some adjustments to the analyses have been made to make them less 
sensitive to background. 

This paper is organized as follows. In Section 2 we describe the data used and the Monte Carlo 
simulation of signal and background processes. In Section 3 we describe the selections for the processes 
ZZ — > , ZZ — > l + £~vV, ZZ — > qqfl>, ZZ — ► qq£ + ^~ and ZZ — > qqqq, where £ + £~ denotes a 




e+ Z e + Z 



Figure 1: NC2 Feynman diagrams for the process e + e — ► ZZ leading to a final state with four 
fermions. 
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Energy point label 
(GeV) 


Mean center-of-mass energy 
(GeV) 


Approximate integrated luminosity 
(pb- 1 ) 


192 


191.59 ±0.04 


29 


196 


195.53 ± 0.04 


73 


200 


199.52 ± 0.04 


74 


202 


201.63 ±0.04 


37 


205 


204.80 ± 0.05 


80 


207 


206.53 ± 0.05 


140 



Table 1: Luminosity- weighted mean center-of-mass energies and typical integrated luminosities for the 
data presented here. The integrated luminosities for each channel vary slightly since they depend on 
the status of different detector elements. The channel dependent luminosities are given in Tables 5 - 
10. 

charged lepton pair of opposite charge and qq any of the five lightest quark-antiquark pairs. We also 
consider the processes ZZ — > bb£ + £~, ZZ — > bbz^F and ZZ — ► qqbb using b-tagging methods similar to 
those developed for the OPAL Higgs search [10]. In Section 4 the selected events are used to measure 
the Z-pair cross section. In this section we also test our measured values of the Z-pair cross section 
and the value of the branching ratio BR(Z — ► bb) measured in Z-pair decays for consistency with the 
Standard Model prediction. At LEP, Z-pair events with b quarks in the final state form the principle 
background to the search for Higgs bosons produced in association with a Z boson. In Section 5 
anomalous neutral-current triple gauge couplings are constrained using the measured cross sections 
and kinematic variables combined in an optimal observables (OO) method. Finally, in Section 6 limits 
are placed on low scale gravity theories with large extra dimensions. 

2 Data samples and Monte Carlo simulation 

The OPAL detector 1 , trigger and data acquisition system used for this study are described in Refer- 
ences [11-15]. 

The approximate integrated luminosities and luminosity-weighted mean center-of-mass energies [16] 
for the data used in this analysis are given in Table 1. The actual integrated luminosities used for each 
final state vary due to differing detector status requirements. These luminosities were measured using 
small-angle Bhabha scattering events recorded in the silicon-tungsten luminometer [15, 17, 18] and the 
theoretical calculation given in Reference [19]. The overall uncertainty on the luminosity measurement 
amounts to less than 0.5% and contributes negligibly to our cross-section measurement error. 

Selection efficiencies and backgrounds were calculated using a simulation [20] of the OPAL detector. 
The simulated events were processed in the same manner as the data. We define the ZZ cross section 
as the contribution to the total four-fermion cross section from the NC2 Z-pair diagrams shown 
in Figure 1. All signal efficiencies given in this paper are with respect to these Z-pair processes. 
Contributions from all other four-fermion final states, including interference with NC2 diagrams, are 
considered as background. For studies of the signal efficiency we have used grc4f [21] and YFSZZ [22] 
with PYTHIA [23] used for the parton shower and hadronization. YFSZZ only generates the NC2 
diagrams. The grc4f Monte Carlo includes four-fermion background processes and the interference of 
background and the NC2 ZZ signal. Weights, based on the grc4f matrix elements, are assigned to each 

J OPAL uses a right-handed coordinate system in which the z axis is along the electron beam direction and the x axis 
points towards the center of the LEP ring. The polar angle, 0, is measured with respect to the z axis and the azimuthal 
angle, 4>, with respect to the x axis. 
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grc4f Monte Carlo event for the event to originate from NC2 ZZ signal, four-fermion background and 
interference between the four-fermion background and the NC2 ZZ signal. In all of the fits performed 
in this analysis we use the ZZTO calculation [24] with the Gf renormalization scheme for integrated 
Standard Model cross sections. To check the simulation of anomalous couplings in the YFSZZ Monte 
Carlo we use EEZZ [25]. 

Backgrounds are simulated using several different generators. KK2f [26] (with PYTHIA used for 
the parton shower and hadronization) is used to simulate two-fermion final states such as e + e~ — > 
Z*(nj) — > qq(«7) and c + c~ — > 7* (717) — > qq(n7), where (717) indicates the generation of one or more 
initial-state radiated photons. HERWIG [27] and PYTHIA are used as checks for these final states. 
These two-fermion samples include gluon radiation from the quarks, which produce qqg, qqqq and 
qqgg final states. KK2f is also used to simulate muon and tau pair events and Bhabhas are modeled 
using BHWIDE [28] and TEEGG [29]. The grc4f generator, with the contribution exclusively due 
to NC2 diagrams removed, is used to simulate other four-fermion background. KORALW [30] and 
EXCALIBUR [31] are used as checks of the four-fermion background. 

Multiperipheral ("two-photon") processes with hadronic final states are simulated by combining 
events from PHOJET [32], for events without electrons 2 scattered into the detector, and HERWIG [27], 
for events with single electrons scattered into the detector. Two-photon events with both the electron 
and positron scattered into the detector, which are only a signficant background for the qqe + e~ 
final state, are simulated with TWOGEN [33]. The Vermaseren [34] generator is used to simulate 
multiperipheral production of the final states e + e~l + £~ . 

To avoid background from four-fermion final states mediated by a Z boson and a virtual 7*, our 
selections were optimized on Monte Carlo samples to select events with candidate Z-boson masses, 
mi and 7712, that satisfy mi + mi > 170 GeV and \m\ — m^\ < 20 GeV. Above 190 GeV more than 
90% of the events produced via the NC2 diagrams are contained in this mass region. Events from the 
NC2 diagrams dominate in this mass region except for final states containing electron pairs. However, 
inside the acceptance of the electromagnetic calorimeters 3 , the backgrounds from two-photon and 
electroweak Compton scattering (e7 — > eZ) processes [35] to Z final states with one Z decaying to 
electron pairs is reduced to less than 5% of the expected Z-pair cross section and the NC2 diagrams 
again dominate. 

3 Event selection 

The OPAL selections cover all ZZ final states except vvvv and t + t~vv. In hadronic final states, the 
energies and directions of the jets are determined using tracks to reconstruct charged particles and 
electromagnetic and hadronic calorimeter clusters to reconstruct neutral particles. The correction for 
unavoidable double counting due to the overlap of calorimeter energy deposited from charged and 
neutral particles is described in Reference [10]. 

In the qqe + e~, qq/j + /J~ and qqqq analyses, four-constraint (4C) and five-constraint (5C) kinematic 
fits are used. The 4C fit imposes energy and momentum conservation. In the 5C fit the added 
constraint requires the masses of the two candidate Z bosons to be equal to one another. For final 
states with either Z boson decaying to a tau pair, the direction of each tau lepton is approximated by 
the reconstructed particles that identify the tau. The energies and total momenta of the tau leptons 
are obtained by leaving the reconstructed direction of the four fermions fixed and scaling the energy 
and momentum of each of the fermions to obtain energy and momentum conservation. The scaled 

2 In this paper reference to a specific fermion, such as an electron, also includes the charge conjugate particle, in this 
case the positron. 

3 The acceptance of the electromagnetic calorimeters for electrons is approximately given by | cos# c | < 0.985, where 
9 e is the polar angle of the electron. 
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values of the tau momentum and energy are then used in the subsequent steps of the analysis. In the 
qqr + r~ final state, subsequent kinematic fits are therefore effectively 2C and 3C fits. 

The selection procedures for events containing charged leptons or neutrinos (£ + £~ £ + t~ ,£ + £~ uV \ 
qq^ + .£~ and qqz/F) are largely unchanged from the analysis used at 183 GeV and 189 GeV in Refer- 
ence [3] , except that the likelihood functions and some of the cuts have been optimized for the higher 
energies. More significant changes have been made to the qqqq and qqbb selections. 

3.1 Selection of ZZ £+£-£+£- events 

Z-pairs decaying to final states with four charged leptons (£ + £~ £ + £~) produce low multiplicity events 
with a clear topological signature that is exploited to maximize the selection efficiency. The £ + £~£ + £~ 
analysis begins by selecting low multiplicity events (less than 13 reconstructed good tracks and less 
than 13 electromagnetic clusters) with visible energy of at least 0.2-^/s and at least one good track 
with momentum of 5 GeV or more. Good tracks are required to be consistent with originating from 
the interaction point and to be composed of space points from at least half of the maximum possible 
number of central tracking detector (Jet Chamber) sense wires. 

Using a cone algorithm, the events are required to have exactly four cones of 15° half angle each 
containing between 1 and 3 tracks. Cones of opposite charge are paired 4 to form Z boson candidates. 

Lepton identification is only used to classify events as background or to reduce the number of cone 
combinations considered by preventing the matching of identified electrons with identified muons. 
Electrons are identified on the basis of energy deposition in the electromagnetic calorimeter, track 
curvature and specific ionization in the tracking chambers. Muons are identified using the association 
between tracks and hits in the hadron calorimeter and muon chambers. 

To reduce background from two-photon events with a single scattered electron detected, we elim- 
inate events with forward-going electrons or backward-going positrons with the cuts cos# c - < 0.85 
and cos# e + > —0.85. Here 9 e - ( 9 e +) is the angle of the electron (positron) with respect to the incom- 
ing electron beam. Background from partially reconstructed qq(n7) events and two-photon events is 
reduced by requiring that most of the energy not be concentrated in a single cone. Defining E V - 1S as 
the total visible energy of the event and -E^ne as ^ ne energy contained in the most energetic cone we 
require E vis - E™% > 0.2 ^s. 

The invariant masses of the lepton pairs are calculated in three different ways which are motivated 
by the possibility of having zero, one or two tau pairs in the event. The events are first classified 
according to the assumed number of tau pairs in the event: 

1. Events without tau pairs can be selected by requiring a high visible energy. Therefore, events 
with E v i s > 0.9-y/i are treated as e + e~e + e~, e + e~ or [i + [i~fi + fi~ events. We also treat all 
events with |cos# m i ss | > 0.98 (6> m i ss is the polar angle associated with the missing momentum 
in the event) as e + e~e + e~, e + e~ / u + / u~ or events to maintain efficiency for Z-pairs 
with initial-state radiation. As there are no missing neutrinos in these events, the mass of each 
cone-pair combination is evaluated using the measured energies and momenta of the leptons. 

2. Events failing (1) with a cone-pair combination that has energy exceeding 0.9mz are tried as 
an e + e _ T + r~ or ^ + /i~r + T~ final state. The mass of the tau-pair system is calculated from the 
recoil mass of the presumed electron or muon pair. 

4 Two-track cones are assigned the charge of the more energetic track if the momentum of one track exceeds that of 
the other by a factor of 4. Events with a cone which fails this requirement are rejected. 
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3. Events failing (1) with a cone-pair combination failing (2) are treated as t + t~t + t~ final states. 
The momenta of the tau leptons are determined using the scaling procedure described in the 
introduction to Section 3. The invariant masses of the cone pairs are evaluated using the scaled 
momenta. 

In any event where the lepton identification allows more than one valid combination, each com- 
bination is tested using invariant mass cuts. For events with one or more combinations of cone- 
pairings satisfying |m z — mu\ < 0.1m z and |m z — mgip\ < 0.1m z the one with the smallest value of 
(m z — mu) 2 + (m z — m^) 2 is selected for further analysis. In the other events, the combination with 
the smallest value of |m z — mu\ or |m z — m^'| is selected. These requirements maintain efficiency 
for signal events with a single off-shell Z boson. The final event sample is then chosen with the re- 
quirement mu + m£'c > 160 GeV and |m« — m^/| < 40 GeV. The signal detection efficiency, averaged 
over all £ + £~£ + £~ final states is approximately 56%. The efficiency for individual final states ranges 
from about 30% for t + t~t + t~ to more than 70% for These efficiencies have almost 

no dependence on center-of-mass energy. In Table 2 (line a) we give the efficiency, background and 
observed number of events. The errors on the efficiency and background include the systematic uncer- 
tainties which are discussed in Section 3.6. The invariant masses of all cone pairs passing one of the 
selections are shown in Figure 2a. A total of four £ + £~£ + £~ events is observed between 190 GeV and 
209 GeV with an expected background of 1.08 ± 0.27 events. The remaining background is dominated 
by four-lepton events which are not from the NC2 diagrams that define Z-pair production. 

3.2 Selection of ZZ — > i+l'iyV events 

The selection of the e + e~vV and ji + n~vv final states is based on the OPAL selection of W pairs 
decaying to leptons [36]. The mass and momentum of the Z boson decaying to vv are calculated using 
the beam-energy constraint and the reconstructed energy and momentum of the Z boson decaying to 
a charged lepton pair. A likelihood selection based on the visible and recoil masses as well as the polar 
angle of the leptons, is then used to separate signal from background. 

The e + e~vV selection starts with OPAL W-pair candidates where both charged leptons are clas- 
sified as electrons. Each event is then divided into two hemispheres using the thrust axis. In each 
hemisphere, the track with the highest momentum is selected as the leading track. The sum of the 
charges of these two tracks is required to be zero. The determination of the visible mass, m v i s , and 
the recoil mass, m Teco n, is based on the energy as measured in the electromagnetic calorimeter and the 
direction of the leading tracks. 

The likelihood selection uses three variables: QcosO, where is the angle of the track with the 
highest momentum and Q is its charge, the normalized sum of visible and recoil masses (m v ; s + 
fyhecmi) / \fs and the difference of visible and recoil masses, m v - m — m reco ii- The performance of the 
likelihood selection is improved with the following preselection: —25 GeV < m v - ls — m reC oii < 15 GeV 
and (m v i s + m re coii) > 170 GeV. Two events with C e + e - vV > 0.60 are selected. See Table 2 (line b) and 
Figure 3a. The expected background is primarily from W pairs and amounts to 1.28 ± 0.15 events. 

The fi + [L~vV selection starts with the OPAL W-pair candidates where both charged leptons are 
classified as muons. The selection procedure is the same as for the e + e~vV final states except that m v i s 
and m rcco ii are calculated from the momentum of the reconstructed tracks of the Z boson decaying to 
muon pairs. The likelihood preselections —25 GeV < m v ; s — m rcco ii < 25 GeV and (m v i s + m reco ii) > 
170 GeV are applied. No event with C^+^- uV > 0.60 is selected while 4.30 ± 0.39 are expected (see 
Table 2 (line c) and Figure 3b). The probability to observe no event when 4.3 ± 0.39 are expected is 
approximately 1.4%. 
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Figure 2: (a) Invariant masses of £ + £~£ + l~ cone pairs. Invariant mass pairs for the (b) qqe + e~ (c) 
qq/i + /^ _ and (d) qqr + T~ data. The dashed lines show the final invariant mass cuts. In contrast to 
Reference [3], the qqe + e~, qq/i + /x~ unconstrained invariant masses are plotted after the requirement 
that the 5C fit mass, with both quark and lepton-pair masses constrained to be the same, exceeds 
85 GeV. The qqr + T~ unconstrained invariant masses are shown after a similar cut on the 3C fit mass. 
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cut used to select ZZ events in each case. 



9 



3.3 Selection of ZZ — > qq£+£ events 



The lepton pairs in the qqe + e~ and qq> + fi~ final states have a distinctive signature making possible 
selections with high efficiencies and low background contaminations. In the qqr + r~ final state, the 
decay of the tau leptons produces events which are more difficult to identify. The identification of 
this final state exploits the missing momentum and missing energy carried away by the neutrinos 
produced in the decay of the tau lepton. In comparison to our previous publication [3], the qqe + e~ 
and qqfi + fi~ selections are almost unchanged, but the analysis of qq~T + T~ events with b-tags has been 
simplified. The qqe + e~ and qq~[i + fi~ selections presented here are also similar to those used in [37] 
but are optimized specifically for Z-pair final states. 

3.3.1 Selection of ZZ — > qqe+e - and ZZ — > qq/x+/x _ events 

The selection of qqe + e~ and qq/x + /x~ final states requires the visible energy of the events to be greater 
than 90 GeV and at least six reconstructed tracks. Among all tracks with momenta greater than 2 GeV, 
the highest momentum track is taken as the first lepton candidate and the second-highest momentum 
track with a charge opposite to the first candidate is taken as the second lepton candidate. Using the 
Durham [38] jet algorithm, the event, including the lepton candidates, is forced into four jets and the 
jet resolution variable that separates the three-jet topology from the four-jet topology, 2/34, is required 
to be greater than 10 -3 . Excluding the electron or muon candidates and their associated calorimeter 
clusters, the rest of the event is forced into two jets. The 4C and 5C fits to the two lepton candidates 
and the two jets are required to converge. 5 

In the qqe + e~ selection no explicit electron identification is used. Electron candidates are selected 
by requiring the sum of the electromagnetic cluster energies E\ + E2 associated to the electrons to 
be greater than 70 GeV and the momentum of the track associated with the most energetic electron 
candidate to exceed 20 GeV. We also reject the event if the angle between either electron candidate 
and any other track is less than 5°. 

In the qq/x + ^~ selection the muons are identified using (i) tracks which match a reconstructed 
segment in the muon chambers, (ii) tracks which are associated to hits in the hadron calorimeter or 
muon chambers [36], or (iii) isolated tracks associated to electromagnetic clusters with reconstructed 
energy less than 2 GeV. No isolation requirement is imposed on events with both muon tracks passing 
(i) or (ii). Events with at least one muon identified with criterion (iii) are accepted if both muon 
candidates in the event have an angle of at least 10° to the nearest track. We require the sum of the 
momenta of the two leptons to be greater than 70 GeV. 

Z-pair events are separated from Z7* background by requiring the fitted mass of the 5C fit to 
be larger than 85 GeV and the invariant masses mu and m qq obtained from the 4C fit to satisfy 
(m« + m,qq) > 170 GeV and |m« — m qq | < 30 GeV. Figure 2b (2c) shows the distribution of m ee (m^) 
and m qq before the cuts on the masses from the 4C fit. 

After all cuts the mean selection efficiency 6 for qqe + e~ signal events is 60% and depends weakly 
on energy. A total of 17 candidate events is found after all cuts at between 190 GeV and 209 GeV, 
which can be compared to the Standard Model expectation (including a small background) of 14.1±0.4. 
In Table 2 (lines d and e) we give the efficiency, background and observed number of events. The 
errors on these efficiencies and backgrounds include the systematic uncertainties (see Section 3.6). 
The largest source of background after all cuts is from the process e + e~ — ► Z7* — > qqe + e~ and from 
two-photon events. 

5 In the context of this paper, convergence is defined as a fit probability greater than 10 -10 after at most 20 iterations. 
6 Small amounts of feedthrough from other ZZ final states, in this case qqr + T~, are counted as signal. 
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In the qq~/U + ^~ selection, the mean selection efficiency is 73% and varies less than 5% with energy. 
A total of 22 events is observed at between 190 GeV and 209 GeV, which can be compared to the 
Standard Model expectation (including a small background) of 15.5 ± 0.4. (See Table 2 lines f and g). 
The background after all cuts is expected to come mainly from e + e~ — ► Z7* — ► qq/z + // _ events. 

3.3.2 Selection of e+e _ — > ZZ — > qqr+r - events 

The qqr + r^ final state is selected from a sample of events with track multiplicity greater than or equal 
to six. Events which pass all the requirements of the qqe + e~ and qq> + ^~ selections are excluded 
from this selection. The tau-lepton candidates are selected using a neural network algorithm which 
is described in detail in Reference [10]. At least two tau candidates are required. The tau candidate 
with the highest neural network output value is taken as the first candidate. The second candidate 
is required to have its charge opposite to the first candidate and the highest neural network output 
value among all remaining candidates. 

Because of the presence of neutrinos in the final state the missing energy, \fs — E V - 1S , is required to 
exceed 15 GeV while the visible energy of the event, E v i s , is required to exceed 90 GeV. In addition, the 
sum of the momenta of the leading tracks from the tau-lepton decays is required to be less than 70 GeV. 
Since the direction of the missing momentum in signal events will tend to be along the direction of 
one of the decaying tau leptons, the angle a Tim i ss between the missing momentum and a tau-lepton 
candidate is required to satisfy 

QV miss ^ 90° for at least one of the two tau candidates. 

The two hadronic jets are selected in the same way as in the e + e~ — > ZZ — > qqe + e~ selection. The 
initial estimate of the energy and the momenta of the tau candidates is found from the sum of the 
momenta of the tracks associated to the tau by the neural network algorithm and all unassociated 
electromagnetic clusters in a cone with a half angle of 10° around the leading track from the tau 
decay. A 2C kinematic fit that imposes energy and momentum conservation (see the introduction to 
Section 3) is required to converge. A 3C kinematic fit, with the additional constraint of the equality 
of the fermion pair masses is also required to converge. 

Using the network output [10] for each tau candidate, a tau lepton probability, V, is calculated 
taking into account the different branching ratios, sensitivities, efficiencies and background levels for 
1-prong and 3-prong tau-lepton decays. In this paper, we combine two probabilities V\ and Vi, to 
form a likelihood ratio using 

C = ^ (1) 

The likelihood ratio associated with probabilities of the two tau candidates is required to satisfy 
C TT > 0.977. In addition, the common mass of the 3C fit is required to exceed 85 GeV. Using the 2C 
fit masses of the tau pair, m TT , and the quark pair, m qq , as obtained from the kinematic fit, we also 
require m qq + m TT > 170 GeV and \m qq — m TT \ < 60 GeV. 

Four candidate events are found in the data. Figure 2d shows the masses of the candidate events 
before the invariant mass cuts. In Table 2 (lines h and i) we give the efficiencies, backgrounds and 
observed number of events. 



3.3.3 Selection of e+e~ -» ZZ -» hh£+£~ 

The selection of the bb£ + £~ events is based on the qq£ + £~ selections with addition of the algorithm 
described in Reference [10] to identify bb final states. The probabilities that each of the hadronic jets 
is a b jet can be combined to form a likelihood function, £bb> according to Equation (1). Because 
the qqe + e~ and qq~/U + /i~ selections are pure, a relatively loose cut of £bb > 0.2 is used to select the 
bbe + e~ and bb/i + /i~ samples. For the selections with electron and muon pairs there are two classes 
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of events since the selected bbl + l~ events are a subset of the qq£ + £~ events. In Table 2 (lines e and 
g) we give the efficiencies of the b-tagged samples with respect to the expected fraction of bb events. 
The efficiencies for samples without b-tags are given with respect to the hadronic decays without bb 
final states. 

In the data with y/s above 190 GeV we find three candidate bbe + e~ events with 2.49±0.19 expected 
and in the bbfi + fi~ selection we find seven events with 2.59 ± 0.15 events expected. The probability 
to observe seven or more events when 2.59 ± 0.15 are expected is approximately 1.7%. The invariant 
masses of these seven events are consistent with mz ■ 

For the bbr + r~ selection the C TT cut of the qqr + r _ selection is loosened and combined with £t>b 
as follows. C TT and £bb are both required to be greater than 0.1. The bbr + r~ probability for the 
event, Cbbrr, is calculated from Equation (1) with C TT and £bb as inputs and required to exceed 
0.95. After the cut on £t>brT> the remaining cuts of the qqr + T~ selection are applied. None of the 
qq~T + T~ events are identified as bbr + r _ candidates which can be compared with the Standard Model 
expectation of 1.41 ± 0.19. One additional candidate event is found in the sample with the relaxed 
cut on C TT with 0.35 ± 0.10 expected. See Table 2 (lines i and j). 

3.4 Selection of ZZ — > qQvV events 

The qqW selection is based on the reconstruction of the Z boson decaying to qq which produces 
slightly boosted back-to-back jets. The selection uses events with a two-jet topology where both jets 
are contained in the detector. The beam energy constraint is used to determine the mass of the Z 
boson decaying to vv. The properties of the qq decay and the inferred mass of the vv decay are then 
used in a likelihood analysis to separate signal from background. 

Two-jet events are selected by dividing each event into two hemispheres using the plane perpen- 
dicular to the thrust axis. The number of tracks in each hemisphere is required to be four or more. 
The polar angles of the energy-momentum vector associated with each hemisphere, #hemii and #hemi2> 
are used to define the quantity cos#h = ^(cos #hcmii — cos#hcmi2)- Contained events are selected by 
requiring |cos#h| < 0.80. The total energy in the forward detectors and in the forward region of 
the electromagnetic calorimeter (|cos#| > 0.95) is required to be less than 3 GeV. W boson decays 
identified by the OPAL W-pair selection are rejected using the likelihood function for e + e~ — > qq^ 
from Reference [36] which includes an optimization for each center-of-mass energy. Only events with 
£ww < 0.5 are retained. 

An important background to our selection is qq(^7) events with photons that escape detection. We 
discriminate against these events by looking for a significant amount of missing transverse momentum, 
p t . In each event, pt can be resolved into two components, pu, perpendicular to the transverse 
component of the thrust axis and ptj, along the transverse component of the thrust axis, pu depends 
primarily on the reconstructed angles of the jets and therefore is more precisely measured than p t j 
which depends on the energy balance of the jets. We approximate pu as pa = \E\> sin^sin^h- 
Here E^ = \/~s/2 is the beam energy, <f> is the acoplanarity calculated from the angle between the 
transverse components of the momentum vectors of the two hemispheres and sin 9^ = y/l — cos 2 6^. 
The resolution on pa, a Pti , was parameterized as a function of thrust and cos 9^ using data taken at 
the Z resonance. We construct the variable 

Rpu = (Pti-Pu)/v P u, (2) 

which is used as input to the likelihood function described below. Here corresponds to the transverse 
momentum carried by a photon with half the beam energy which just misses the inner edge of our 
acceptance ( p^ = sin(32 mrad)/2). The likelihood function also uses the related variable cos# m i ss , 
the direction of the missing momentum in the event, to discriminate against the qq(^7) events. 

In the final selection of events, we use a likelihood function based on the following five variables: 
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1. the normalized sum of visible and recoil masses (m v i s + m TCCO u) / ^/s , 

2. the difference of visible and recoil masses (m v i s — m rccoi i), 

3. log(y23)> where 2/23 is the jet resolution parameter that separates the two-jet topology from the 
three-jet topology as calculated from the Durham jet algorithm, 

4. cos0 m i ss and 

5. R m . 

The mass variables are useful for reducing background from W-pair production and single W (Weis) 
final states. The jet resolution parameter is useful in reducing the remaining qq£v final states. To 
improve the performance of the likelihood analysis we use only events with |m v i s — m reC oii| < 50GeV, 
(m v i s + m rcco ii) > 170 GeV and R Pti > 1.2. Events are then selected using £ q q U v > 0.5, where C q q V v 
is the likelihood function for the qqW selection. The likelihood distribution of data and Monte Carlo 
simulation is shown in Figure 3c. For the bbvV selection we require, in addition, the b-tag variable of 
Reference [10] to be greater than 0.65. 

The mean efficiency for the qquV selection alone is 32% and does not have a strong dependence on 
s/s. The efficiency includes corrections to the Monte Carlo simulation for the effects of backgrounds 
in the foward detectors and for imperfect detector modeling. The modeling systematic uncertainties 
are discussed below in Section 3.6. A total of 60 candidates is observed in the data to be compared 
with the Standard Model expectation of 64.9 ± 2.9. The efficiencies after considering the results of the 
b-tagging, as well as the number of events selected are given in Table 2 (lines k and 1) . 

3.5 Selection of ZZ — > qqqq events 

As in the previous analysis [3] the hadronic selection consists of a preselection followed by two different 
likelihood selections, one of them aiming at an inclusive selection of fully hadronic Z-pair decays, and 
the other optimized for selecting final states containing at least one pair of b quarks. The main changes 
with respect to the previous analysis are that cuts in the preselection have been modified slightly and 
that the quantities used in the qqbb likelihood have changed. 

3.5.1 Preselection 

The preselection starts from the inclusive multihadron selection described in Reference [17]. The 
radiative process e + e~ — > Z7 — > qq~7 is suppressed by requiring the effective center-of-mass energy 
after initial-state radiation, ^/s 7 , to be larger than 160 GeV. Vs' is obtained from a kinematic fit [17] 
that allows for one or more radiative photons in the detector or along the beam pipe. The final-state 
particles are then grouped into jets using the Durham algorithm [38]. A four-jet sample is formed by 
requiring the jet resolution parameter 2/34 to be at least 0.003 and each jet to contain at least two 
tracks. In order to suppress Z*/7* — ► qq background, the event shape parameter C par [39], which is 
large for spherical events, is required to be greater than 0.27. A 4C kinematic fit using energy and 
momentum conservation is required to converge. A 5C kinematic fit which forces the two jet pairs to 
have the same mass is applied in turn to all three possible combinations of the four jets. This fit is 
required to converge for at least one combination. 

The most probable pairing of the jets is determined using a likelihood discriminant which is based 
on the difference between the two di-jet masses calculated from the results of a 4C kinematic fit, the 
di-jet mass obtained from a 5C kinematic fit, and the \ 2 probability of the 5C fit. 
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3.5.2 Likelihood for the inclusive ZZ — >• qqqq event selection 

We use a likelihood selection with six input variables for the selection of ZZ — > qqqq events. The 
likelihood selection is optimized separately for each energy point. The first variable is the output 
value of the jet pairing likelihood described above. The second variable is determined by excluding 
the jet pairing with the largest difference between the two di-jet masses as obtained from the 4C fit, 
and then evaluating the 5C-fit masses of the remaining two pairings. The one which is closest to the 
W mass is selected and the difference between this 5C-fit mass and the W mass is used in order to 
discriminate against hadronic W-pair events. The third variable suppresses background events with 
radiated photons that shift the value of the mass obtained from the kinematic fit. We use the value of 
v / s' from the preselection. The fourth variable is used to discriminate against Z* /j* events. We use 
the difference between the largest and smallest jet energies after the 4C fit. The final two variables 
are calculated from the momenta of the four jets. They are the effective matrix elements for the QCD 
processes Z*/7* — > qqgg and Z*/7* — > qqqq as defined in Reference [40], and the matrix element for 
the process WW — > qqqq from Reference [31]. 

The cut on the likelihood function has been chosen in order to minimize the total expected relative 
error when including a 10% relative systematic uncertainty on the background rate in addition to the 
expected statistical error. The qqqq likelihood function for the data is shown in Figure 4a. A total of 
206 events is observed which can be compared to the total Standard Model expectation of 201.9 ± 13.6 
which includes a large background of 135.6 ± 13.1 events. In Table 2 (lines m and o) we give the 
efficiency, background and observed number of events. 

3.5.3 Likelihood function for ZZ — ► qqbb event selection 

Jets originating from b quarks are selected using the b-tagging algorithm [10]. We evaluate the 
probability for each of the four jets to originate from a primary b quark. The two highest probabilities 
are then used as input variables for a likelihood function to select ZZ — > qqbb events. In addition 
we use the parameters 2/34, C par , the track multiplicity of the event and the output of the jet pairing 
likelihood function. We also use the fit probability of a 6C fit which forces both masses to be equal to 
the W mass. Finally we use the variable | cos 6>bb — cos qq |, where #bb is the opening angle between the 
most likely b jets, and 6> qq is the opening angle between the remaining two jets 7 . The qqbb likelihood 
function for the combined data taken at y/s > 190 GeV is shown in Figure 4b. A total of 45 events is 
observed which can be compared to the total Standard Model expectation of 43.3 ± 1.8. In Table 2 
(lines n and o) we give the efficiency, background and observed number of events. 

3.6 Systematic uncertainties 

Systematic uncertainties have only a modest effect on our final result because of the large statistical 
error associated with the small number of Z-pair events produced at LEP 2. 

Detector systematic uncertainties for the qq£ + £~ and £ + £~ £ + £~ selections without r-pairs in the 
final state are small because of the good separation of signal and background. In the qq£ + £~ channels 
possible mismodeling of the detector resolution is accounted for by smearing polar and azimuthal jet 
angles with a Gaussian width of 1° and the energies by 5%. This gives systematic shifts smaller than 
5%. These shifts are included in the systematic uncertainties on the efficiencies. In the qq~T + T~ selec- 
tion, the additional systematic uncertainties in the efficiencies are determined by overlaying hadronic 
and tau decays taken from Z resonance data giving the dominant contribution to the total systematic 

7 In contrast to our previous publication [3] , the probability of the 5C mass fit that constrains one of the candidate Z 
bosons to the Z mass, as used in the OPAL Higgs [10] analysis, is not an input to the likelihood function. 
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Figure 4: The combined distributions of likelihood discriminants used for all energies above 190 GeV. 
(a) e + e~ — ► ZZ — ► qqqq selection and (b) e + e~ — ► ZZ — > qqbb selection. The dotted line and arrow 
show the position of the likelihood cut used to select ZZ events in each case. 
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Table 2: Observed number of events, n Q bs, the total Standard Model expectation, usm, the expected 
number of Z-pairs, nzz, background expectation, nb a ck> an d efficiencies, e c h a n> f° r the combined data 
sample. Bzz is the product branching ratio for the final state which is calculated directly from Z 
resonance data [6]. An overbar is used to indicate that events from a particular selection are rejected. 
Note that the efficiencies for selections with b-tags are given relative to the fraction of hadronic final 
states which contain a Z boson decaying to bb. For selections of events with hadronic final states, 
but without b-tags, the efficiencies are relative to those expected hadronic final states which do not 
include a Z boson decaying to bb. The errors on all quantities include contributions from the systematic 
uncertainties described in Section 3.6. 

uncertainty of 6.2%. In the £ + £~£ + £~ final state the largest effect (3%) is from the modeling of the 
multiplicity requirement which is important for final states containing r-pairs. 

In the final states with neutrinos, £ + £~vV and qqVTy, tight cuts are needed to separate signal and 
background. In these selections detector effects can best be studied by comparing calibration data 
taken at the Z resonance with a simulation of the same process. In these cases, we add additional 
smearing to the total energy and momentum of the simulated events to match data and simulation on 
the Z-resonance. We then apply the same smearing to the signal and background Monte Carlo simu- 
lations for the Z-pair analyses. The reported efficiencies and backgrounds are accordingly corrected. 
The full difference is used as the systematic uncertainty in these cases. These differences give relative 
systematic uncertainties on the efficiency of 2.5% for the e + e~uV final state, 5.2% for the ^ ' [i~vv 
final state and 3.8% for qqV77 final state. 

In the qqqq inclusive analysis the sensitivity to the detector description of jet eneriges is much 
reduced by the use of kinematic fits. In this case the systematic uncertainties are determined by 
smearing the jet energies by an additional 5% and the jet directions by 1° leading to a relative 
detector systematic uncertainty of 6%. 

Another important detector effect comes from the simulation of the variables used by the OPAL 
b-tag which is discussed in Reference [10]. We allow for a common 5% uncertainty on the efficiency 
of the b-tag, consistent with our studies on Z resonance data and Monte Carlo simulation. 

In each channel the signal and background Monte Carlo generators have been compared against 
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alternative generators. In almost all cases the observed differences are consistent within the finite 
Monte Carlo statistics and the systematic uncertainty has been assigned accordingly. 

The largest contribution to the systematic uncertainty from the model dependence of the back- 
ground prediction is for the qqqq and qqbb final states. This uncertainty has been estimated by 
comparing the predictions of KORALW and grc4f for the W-pair background, and the prediction of 
PYTHIA and KK2f for the background from hadronic 2-fermion production. The resulting uncer- 
tainty is 10% on the background for the inclusive selection and 20% for the background of the qqbb 
selection [10]. The 10% uncertainty has been taken to be fully correlated between energy points and 
between the qqqq and qqbb selections. The correlated error on the qqbb selection from uncertainties in 
the background due to b-tagging, has been absorbed into the overall b-tagging efficiency uncertainty. 

In the £ + £~£ + £~ channel, which has a large background from two-photon events, we have compared 
the number of selected events at an early stage of the analysis with the Monte Carlo simulation and 
based our background systematic uncertainty on the level of agreement. This results in 20% systematic 
uncertainty on the background. 

The efficiencies and backgrounds given for each energy point (see the Appendix and the summary 
Table 2) include the systematic uncertainties described in this section as well as errors from finite 
Monte Carlo statistics. Given the large statistical error from the small Z-pair cross section we have 
employed the following conservative scheme to account for possible correlations among the systematic 
uncertainties: 

1. The common correlated error from detector systematics, hadronization and kinematic distri- 
butions of the Z-pair events is taken to be 3%. This is based primarily on the studies of the 
smearing of reconstructed jet and lepton directions. 

2. The common systematic uncertainty for any channel including a b-tag is taken to be 5%. 

3. As described above an uncertainty of 10% of the backgrounds is taken to be fully correlated 
among all qqqq and qqbb channels. 

4. Based on generator level comparisons of kinematic distributions we conservatively assume a 2% 
correlated systematic uncertainty for the modeling of any changes to the physics description of 
the Z process when new physics is switched on. This is in addition to an overall uncertainty of 
2% for the Standard Model total cross section prediction of the ZZTO calculation. 

4 Cross section and branching ratios 

At each center-of-mass energy point, information from all of the analyses is combined using a maximum 
likelihood fit to determine the production cross section for e + e~ — > ZZ. The information which was 
used in the fits, as well as the Standard Model prediction for Z-pair production, is summarized in 
Table 2. The details for each energy point are given in the Appendix. For each channel the tables give 
the number of events observed, n bs, the Standard Model prediction for all events, tism, the expected 
signal, razz, the expected background, raback, the efficiency e c h an , and the integrated luminosity, Lint- 
Bzz is the branching ratio of Z-pairs to the given final state, calculated from Z resonance data [6]. 
In the table we give the overlap between the b-tag and non-b-tag analyses. Possible overlap between 
qqqq and qq£ + £~ has been studied, and found to be an order of magnitude smaller than the overlap 
of qqqq and qqbb and has therefore been ignored. 

The cross section at each energy is determined with a maximum likelihood fit using Poisson proba- 
bility densities convolved with Gaussians to describe the uncertainties on efficiencies and backgrounds. 
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Figure 5: The OPAL measurements of the NC2 Z-pair production cross section. The middle of the 
shaded band shows the Standard Model prediction of the ZZTO calculation. The band indicates 
a theoretical uncertainty of ±2%. The error bars on the data show the combined statistical and 
systematic uncertainty, the horizontal marks indicate the extent of the statistical errors. Data from 
y/s near 183 GeV and 189 GeV are from Reference [3]. 
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Figure 6: The OPAL measurements of the NC2 Z-pair differential cross section d<7zz/d| cos#z|- The 
curves show the prediction of the YFSZZ Monte Carlo for the differential cross section, which has 
been normalized to agree with the total cross section of the ZZTO calculation. The error bars show 
the combined statistical and systematic uncertainty. The arrow shows the 95% confidence level upper 
limit for the one case where the minimum of the negative log likelihood function was at zero. 
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The expected number of events in each channel, fx e , as function of the Z-pair cross section, azz, is 
given by 

(3) 

The efficiencies, e c ham include the effects of off-shell Z bosons that are produced outside of our kine- 
matic acceptance. The correlated systematic uncertainties described in the last section are imple- 
mented by introducing additional parameters which are constrained with Gaussian probability densi- 
ties given by the size of the systematic uncertainty. Our main results, the NC2 Z-pair cross sections 
obtained from the fits, are 
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The first error is statistical and the second error is systematic. The systematic errors are obtained from 
the quadrature difference of the errors when running the fit with and without systematic uncertainties. 
The comparison of these measurements with the ZZTO prediction is shown in Figure 5. The results 
are consistent with the Standard Model prediction. 

We have also analyzed our data in four bins of | cos 0z\ where Qz is the polar angle of the Z bosons. 
Except for channels with one Z boson decaying to neutrinos, the value of | cos#z| is determined from 
kinematic fits which assume no initial-state radiation. In the t^i^vv and qqVF channels, the direction 
is determined from the reconstructed direction of the visible Z. The comparison of the expected and 
observed differential cross sections is shown in Figure 6. The data in Figure 6 have been corrected for 
small amounts of bin migration assuming the Standard Model prediction. 

In order to check the consistency of our result with the Standard Model, we have performed a 
maximum likelihood fit in which the Standard Model ZZTO prediction is scaled by an overall factor 
R. The fit is based on individual measurements for each channel given in the Appendix in Tables 5 
to 10 as well as the data presented in Reference [3]. The treatment of the correlated systematic 
uncertainties is outlined in Section 3.6. The fit yields 

where the error includes the correlated experimental systematic uncertainty (3%) and the theoretical 
uncertainty on the ZZTO prediction (2%). In a separate fit, we allowed the branching ratio of Z 
bosons to b quarks, BR(Z — ► bb), to be a free parameter in the fits at each of the energy points. The 
results are shown in Figure 7. The BR(Z — ► bb) values tend to lie above the measured value from 
LEP 1 data; combining all energies, the average value is 0.196±0.032, which is 1.3 standard deviations 
above the LEP 1 measurement of 0.1514±0.0005 [6]. 

5 Limits on anomalous triple gauge couplings 

Limits on anomalous triple gauge couplings were set using the measured cross sections and kinematic 
information from an optimal observable (OO) method for the qqqq, qqbb and qc[l + £~ selections. In 
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Figure 7: Determination of Br(Z — ► bb) from ZZ events. The line shows the measured value from 
LEP 1 data. Only statistical errors are shown. The data have a common systematic uncertainty 
which amounts to approximately 5% of the expected branching ratio. The 189 GeV point is from 
Reference [3]; at 183 GeV the ZZ event samples were too small to allow a meaningful measurement to 
be obtained. 
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this study we vary the real part of the ZZZ and ZZ7 anomalous couplings parameterized by /<f zz , 
/ zzz , / Z 27 and /^ Z7 as defined in Reference [7] and implemented in the YFSZZ Monte Carlo. In most 
cases, the real parts of each coupling were varied separately with all others fixed to zero. 

The 00 analysis used here is described in detail in Reference [41]. Since the effective Lagrangian 
used to describe the anomalous couplings is linear in the couplings, the resulting differential cross 
section is parabolic and can be parameterized, for a single non-zero coupling aj, as 

^ = s^(n) + a 4 1 \n) + a^ 2 \n) (4) 

where ft is the phase-space point based on the four-momenta of the four out-going fermions from the 
Z decays. The optimal observables are defined as 

o\ ee sf\n)/s (0 \n). { ' 



To reduce the dependence of our analysis on the tails of the 0\ and 0\ distribution, extreme values 
of 0\ and 0\ were rejected when calculating the mean values. These cuts were chosen based on the 
expected distribution of 0\ and 0\ calculated from Monte Carlo and typically reject a few percent of 
the expected events. 

The expected average values of the first order optimal observable, (0\), and the second order 
optimal observable, (O2), were calculated using the YFSZZ matrix element to reweight accepted 
events from the signal Monte Carlo. The same cuts on extreme values of 0\ and O l 2 were used in data 
and Monte Carlo simulation. The parameterization used was of the form 

Po + p\a +P2a 2 
do + d\oi + (L2OL 1 

(6) 

do + d\a + c?2« 2 

where a is the value of the anomalous coupling and the coefficients are determined separately for 
each selection, energy and type of anomalous coupling. Note that the mean values of the optimal 
observables are normalized to the observed number of events and do not depend on the cross section 
which is considered separately below. This can be seen in Equation (6) by noting that the polynomial 
do + d\a + c?2« 2 parameterizes the change in cross section with anomalous coupling. The average 
value of the optimal observable for the exclusive qqe + e~, qq[i + fi~ and three categories of qqqq events 
( qqqq & qqbb, qqbb & qqqq and qqbb & qqqq ) were calculated from the data at each energy between 
190 GeV and 207 GeV. The statistical uncertainties on these values were parameterized with a co- 
variance matrix determined from high statistics Standard Model signal and background Monte Carlo 
simulations. The covariance matrix was then scaled to match the number of events observed in the 
data. 

The systematic uncertainties were taken into account by varying the modeling of the Standard 
Model signal and the backgrounds. Any deviations were interpreted as a systematic error and included 
in the covariance matrix. The error associated with the physics simulation of the signal was determined 
by comparing the values of {0\) and {O2) determined from reweighting the YFSZZ and grc4f Monte 
Carlo event samples. The error in the background determination was evaluated by using the alternative 
background simulations described in Section 3.6. These differences were taken to be fully correlated 
among energies for a given final state. The systematic uncertainties associated with the accuracy of 
the event reconstruction were evaluated by applying additional smearing to the energies and angles of 
the reconstructed jets and leptons. The resulting contributions to the covariance matrix were taken 
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Coupling 


95% C.L. lower limit 


95% C.L. upper limit 


fZZZ 
J A 


-0.45 


0.58 


fZZZ 
J5 


-0.94 


0.25 


J.ZZ7 
J A 


-0.32 


0.33 


f ZZ7 

h 


-0.71 


0.59 



Table 3: The 95% confidence level limits on possible anomalous triple gauge couplings. These results 
include all channels and all energies, including the OPAL data from 183 GeV and 189 GeV presented 
in Reference [3]. 

to be fully correlated among all five final state selections. The resulting x 2 functions, converted to 
likelihood curves, are shown as dashed lines in Figure 8. 

For those channels used in the OO analysis, we use in addition only the cross section integrated over 
I cos 6>z|- For channels that are not used in the OO analysis, the cross section in four bins of | cos 9%\ is 
used. The ZZTO calculation is used for the prediction of the Standard Model integrated cross section. 
The change in the Z-pair cross section as a function of | cos 6z | and the values of the anomalous 
couplings is parameterized using the YFSZZ Monte Carlo. In addition, the selection efficiencies for all 
final states are parameterized as a function of the couplings and for each bin in | cos 6z | • An uncertainty 
of 10%, dominated by Monte Carlo statistical errors, is assigned to the correction we applied to all 
of these efficiencies. The resulting likelihood curves based on cross-section measurements are shown 
as dotted lines in Figure 8. These cross-section results include data at 183 GeV and 189 GeV from 
Reference [3]. 

Combining the likelihood associated with the x 2 fit and the likelihood curve from the cross-section 
fit, the 95% confidence level (C.L.) limits on the anomalous couplings were obtained and are given in 
Table 3. The combined likelihood curves are shown as solid lines in Figure 8. 

Equation (4) can easily be extended to the case of two non-zero couplings and the constraints can 
be derived for pairs of couplings [41]. Note that for the two-dimensional OO fits, the event sample 
is slightly different from that used in the one-dimensional OO fits, as cuts are placed on the values 
of optimal observables for more than one coupling. The constraint from the expected cross sections 
can also be calculated for two non-zero couplings. The resulting 95% confidence level contours in the 
/4 and /s plane, including both OO and cross-section constraints, are shown in Figure 9. The 95% 
confidence level corresponds to a change in log likelihood of 3.0 from the minimum value. 



6 Limits on low scale gravity theories with large extra dimensions 

We have also examined the possible effects of low scale gravity (LSG) theories with large extra di- 
mensions [42] on the Z-pair cross section. In the LSG theories considered here, gravity is allowed to 
propagate in D = 4 + n dimensions, while all other particles are confined to four dimensional space. 
Newtonian gravity in three spatial dimensions holds if 

M| lanck oc R n M% +2) (7) 

where Mpi anc k is the Planck scale in the usual four dimensions, Mp is the Planck scale in the D- 
dimensional space, and R is the compactification radius of the n extra dimensions. The case of n = 1 
is excluded by cosmological observations. For the case of n = 2, severe constraints are imposed by 
studies of the gravitational interaction on sub-millimeter distance scales [43,44]. 
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Figure 8: The negative log likelihood as a function of the real part of the four anomalous couplings: 
dotted, constraint from cross-section information; dashed, constraint from optimal observables analy- 
sis; and solid, sum. The 95% confidence level corresponds to a change in log likelihood of 1.92 from 
the minimum value. These results are dominated by the 190 — 209 GeV data, but also include the 
OPAL data from 183 GeV and 189 GeV presented in Reference [3]. 
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Figure 9: The 95% confidence level contours, corresponding to a change of 3.0 in the log likelihood 
from the minimum, for (a) the / ZZ7 — ff zz plane and (b) the / ZZ7 — /f zz plane. These results 
include all channels and all energies, including OPAL data from 183 GeV and 189 GeV presented in 
Reference [3] . The star shows the location of the minimum, the solid dot the Standard Model. 



Because gravity can propagate in the extra dimensions, the amplitude for Z-pair production has 
leading order contributions from the s-channel exchange of Kaluza-Klein graviton states. The effec- 
tive theory at LEP 2 energies is sensitive to the mass scale M s that is used to regulate ultraviolet 
divergences. M s is expected to be of the same order of magnitude as Mp. Using this mass scale, 
M s , a contribution to the Born level amplitude that is proportional to A/Mg is obtained. Here A is 
an effective coupling. In this paper we use the definition of M s given in Reference [45] for our main 
result. 

In order to search for deviations from the Standard Model compatible with LSG theories, the 
differential Z-pair cross section is written as 

tSt = ^; £ ^ z ( s ,cos0) C^{8)M^{K,e + ,e-,a,t) + C^(a)M g ^{K,e + ,e- 1 8 l t) 

Q COS t/2 7*S t- c . c 

(8) 

where a is the fine structure constant, s and t are the Mandelstam variables, k is the electron helic- 
ity, e + and e_ are the polarizations of the Z boson, (3 is the Z velocity, and M^. n (K,e + , e_, s,t) and 
-^gravity (k, e + , e_ , s, t) are the Born level SM and gravity matrix elements [9] . Here -M grav it y (k, e + , e_ , s, t) 
is proportional to -pj. The factor Cr z (s, cos 6) corrects for the effects of the finite width of the Z. The 

factors Cj^(s) and C^g^s) correct for initial-state radiation. 
The factors C§^(s), and Cg^s) are given by 



T YFSZZ 



C ~ \ ^YFSZZ ( 9 ) 
V °Born 



where <7isr SZZ is the YFSZZ cross section with ISR and o"Born ZZ ^ s ^ ne Born level YFSZZ cross section. 
The ISR correction for Standard Model Z production, Cfgj^s), is estimated by running YFSZZ with all 
anomalous couplings off. For the case of gravity, Cfg^ v (s) is estimated by using a choice of anomalous 
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couplings which causes the s-channel to dominate. These factors range from 0.84 at 183 GeV to 0.93 
at 207 GeV. Cf s ^(s) and Cfgg(s) differ from each other by no more than 5%. 

The effect of the finite Z width depends on | cos#z| and is obtained from 



Cr z (s, cos#) = 



^YFSZZ 
°|cos0 z | 

\ a \cose z \ 



(10) 



where a Fz is the prediction of Equation (8) with Cr z (s,cos#) = 1 and A = 0. At center-of-mass 
energies above 195 GeV, Cr z (s,cos0) is within 5% of unity. 

In our fit, Equation (8) is used only to compute the expected difference in cross section between the 
Standard Model with and without LSG switched on. As in the case of the anomalous couplings, the 
prediction of ZZTO for the Standard Model cross section is used and the expected angular dependence 
is taken from YFSZZ. 

A maximum likelihood fit is performed in four bins of | cos6>| for each selection and energy point. 
We also use previously published [3] data binned in four bins of | cos 9\ from 189 GeV and the integrated 
cross section from 183 GeV. The fit yields 

A/M 4 = 2.6 ±2.3 TeV~ 4 . 

The likelihood function is shown in Figure 10 and is approximately parabolic. Note that with the sign 
convention of Reference [45] positive values of A/M 4 give a positive contribution to the cross section. 
Because our measured values of the cross sections are on average slightly larger than the Standard 
Model prediction (see Figure 6 and the value for R given in Section 4), the central value for A/M 4 is 
also positive. 

If we assume that the a priori probability for a theory to be true is uniform in the variable A/M 4 
we can obtain limits on M s for a variety of theoretical approaches using the likelihood curve for 
A/M 4 . First, using the notation of Reference [45], we consider separately theories with A = +1, which 
correspond to positive values of A/M 4 , and with A = —1 which correspond to negative values of A/M 4 . 
In the first (second) case our prior assumes all values of A/M 4 > 0.0 ( A/M 4 < 0.0) are equally likely. 
We also use the approximation that the log-likelihood curve is parabolic. The resulting limits on M s 
are given in Table 4. 

In the approach of Reference [46] a summation of higher order terms is used which eliminates 
the dependence of the prediction on A. Comparing References [45] and [46] limits on M s can be 
determined, as a function of the number of large extra dimensions, n, using the substitution 

A s?" 1 (M s \ . , 

"*-^i+lM "7= ( n ) 



Mr M s " +2 Vv^ 

where the integrals In(^f) are given in Equation B.8 of Reference [46]. As the integrals I n are negative 
for all values of n, limits on M s correspond to the A = +1 case shown in Figure 10. These limits on 
M s as a function of n are given in Table 4. 
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Parameter 


95% C.L. lower limit on M s 


Coupling 


Method of Reference [45] 


A = +l 


0.62 TeV 


A = -1 


0.76 TeV 


Number of extra dimensions 


Method of Reference [46] 


n = 2 


0.92 TeV 


n = 3 


0.82 TeV 


n = 4 


0.73 TeV 


n = 5 


0.67 TeV 


n = 6 


0.62 TeV 


n = 7 


0.59 TeV 



Table 4: The 95% confidence level on M s using the approach of References [45] and [46]. These results 
include all channels and all energies, including the OPAL data from 183 GeV and 189 GeV presented 
in Reference [3]. 




Figure 10: The negative log likelihood as a function of the X/Mf. Limits on M s are determined 
separately for theories with A = +1, which correspond to positive values of X/Mf and A = —1 which 
correspond to negative values of X/Mf. The dashed lines indicate the allowed 95% confidence level 
regions obtained in the two cases. These results are dominated by the 190 — 209 GeV data, but also 
include the OPAL data from 183 GeV and 189 GeV presented in Reference [3] . 



27 



7 Conclusion 



The process e + e~ — > ZZ has been studied at center-of-mass energies between 190 GeV and 209 GeV 
using the final states £ + £~£ + £~, £ + £~vV, qq£ + £~, qqisV, and qqqq. The number of observed events, 
the background expectation from Monte Carlo and the calculated efficiencies have been combined to 
measure the NC2 cross section for the process e + e~ — > ZZ. The measured cross sections for the six 
energy points are: 
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The measurements at all energies are consistent with the Standard Model expectations. Using infor- 
mation from the cross-section measurements and from the optimal observables, no evidence is found 
for anomalous neutral-current triple gauge couplings. The 95% confidence level limits are listed in 
Table 3. We have also derived limits on low scale gravity theories with large extra dimensions which 
are summarized in Table 4. 
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Appendix Selected events at individual energies 

The information which was used in the cross-section fit is summarized in Tables 5 to 10. For each 
channel, the tables give the number of events observed, n Q b s , the Standard Model prediction for all 
events, hsm, the expected signal, nzz, the expected background, nb ac k> the efficiency e c ham an d the 
integrated luminosity, L mt . Note that tism as it appears in the tables is not always exactly n-back + ^zz 
due to rounding. 
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0.28 ±0.03 


0.25 ±0 


02 


0.02 ±0.01 





31 ±0 


03 





037 


28 


6 


i 


qqr + T~ & bb-r + T~ 





0.07 ±0.01 


0.05 ±0 


01 


0.02 ±0.01 





23 ±0 


04 





010 


28 


6 


j 


bbr + r _ & qqr+r - 





0.03 ±0.01 


0.02 ±0 


01 


0.01 ±0.01 





10 ±0 


03 





010 


28 


6 


k 


qqVF & hhvv 


3 


3.06 ±0.22 


1.67 ±0 


14 


1.40 ±0.17 





34 ±0 


03 





219 


28 


6 


1 


qq vv & hhvv 


1 


0.43 ±0.05 


0.36 ±0 


03 


0.07 ±0.03 





27 ±0 


03 





061 


28 


6 


m 


qqqq & qqbb 


27 


17.6 ± 1.4 


3.8 ±0 


3 


13.8 ± 1.4 





56 ±0 


04 





300 


29 


1 


n 


qqbb & qqqq 


2 


0.93 ±0.11 


0.52 ±0 


04 


0.42 ±0.10 





12 ±0 


01 





189 


29 


1 





qqbb & qqqq 


1 


1.68 ±0.13 


1.21 ±0 


07 


0.47 ±0.11 





28 ±0 


02 





189 


29 


1 



Table 5: The 192 GeV data. The description of the entries is given in the caption of Table 2. 
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y/s = 196 GeV 





Selection 


^obs 






^back 


^chan 




Bzz 


Lint 


a 


f+frf+fr 


1 


0.58 ±0.10 


0.37 ±0.01 


0.21 ±0.10 


0.56 ±0 


02 


0.010 


72.5 


b 


c^cTvv 


1 


0.66 ±0.08 


0.37 ±0.05 


0.28 ±0.06 


0.42 ±0 


06 


0.013 


72.7 


c 







0.72 ±0.10 


0.37 ±0.05 


0.36 ±0.08 


0.41 ±0 


06 


0.013 


72.7 


d 


qqe + e~ &bbe+e~ 


2 


2.05 ±0.16 


1.58 ±0.07 


0.47 ±0.14 


0.67 ±0 


03 


0.037 


71.3 


e 


qqe + e~ &bbe + e~ 





0.39 ±0.04 


0.34 ±0.03 


0.05 ±0.02 


0.52 ±0 


05 


0.010 


71.3 


f 


qq/i + /i~ & bb/z+/z _ 


3 


1.97 ±0.07 


1.83 ±0.06 


0.14 ±0.04 


0.77 ±0 


02 


0.037 


71.3 


g 


qq/i + /i _ & bb\x + \x~ 


2 


0.43 ±0.03 


0.42 ±0.02 


0.01 ±0.01 


0.63 ±0 


05 


0.010 


71.3 


h 


qqr + T~ & bbr+r - 





0.81 ±0.07 


0.68 ±0.06 


0.13 ±0.04 


0.29 ±0 


03 


0.037 


71.3 


i 


qqr+r - & bbr+r - 





0.19 ±0.03 


0.17 ±0.03 


0.02 ±0.02 


0.27 ±0 


04 


0.010 


71.3 


j 


bbr + T~ & qqr+T~ 





0.06 ±0.01 


0.05 ±0.01 


0.01 ±0.01 


0.08 ±0 


02 


0.010 


71.3 


k 


qqVF & bbvV 


11 


8.62 ±0.64 


4.79 ±0.43 


3.83 ±0.47 


0.34 ±0 


03 


0.219 


70.3 


1 


qqvV & bbvV 





1.29 ±0.15 


1.00 ±0.10 


0.30 ±0.11 


0.26 ±0 


03 


0.061 


70.3 


m 


qqqq & qqbb 


38 


42.1 ±3.3 


10.3 ±0.7 


31.8 ±3.3 


0.50 ±0 


03 


0.300 


75.4 


n 


qqbb & qqqq 


4 


3.13 ±0.33 


1.89 ±0.14 


1.24 ±0.30 


0.15 ±0 


01 


0.189 


75.4 





qqbb & qqqq 


8 


4.69 ±0.39 


3.31 ±0.21 


1.38 ±0.32 


0.26 ±0 


02 


0.189 


75.4 



Table 6: The 196 GeV data. The description of the entries is given in the caption of Table 2. 



y/s = 200 GeV 





Selection 


T^obs 




nzz 


^back 


^chan 




Bzz 




a 


£+£-£+£- 





0.47 ±0.09 


0.43 ±0.01 


0.05 ±0.08 


0.56 ±0 


02 


0.010 


75.1 


b 


c^cTvv 





0.63 ±0.08 


0.43 ±0.06 


0.21 ±0.06 


0.43 ±0 


06 


0.013 


74.3 


c 


H + H~ vT> 





0.68 ±0.09 


0.39 ±0.06 


0.28 ±0.07 


0.40 ±0 


06 


0.013 


74.3 


d 


qqe + e~ &bbe+e~ 


3 


1.98 ±0.12 


1.69 ±0.08 


0.30 ±0.09 


0.63 ±0 


03 


0.037 


73.7 


e 


qqe + e~ &bbe + e~ 





0.46 ±0.04 


0.41 ±0.03 


0.04 ±0.03 


0.56 ±0 


05 


0.010 


73.7 


f 


qq/i + yU~ & bb,u+/!~ 


5 


2.19 ±0.08 


2.01 ±0.06 


0.18 ±0.05 


0.75 ±0 


02 


0.037 


73.7 


g 


qq/i + /i~ & bb/i + /i~ 


3 


0.43 ±0.03 


0.41 ±0.03 


0.02 ±0.02 


0.56 ±0 


05 


0.010 


73.7 


h 


qqT + T~ & bbr+T~ 


1 


0.90 ±0.08 


0.76 ±0.07 


0.15 ±0.05 


0.28 ±0 


03 


0.037 


73.7 


i 


qqr + T~ & bbr + T~ 





0.23 ±0.03 


0.21 ±0.03 


0.02 ±0.02 


0.28 ±0 


04 


0.010 


73.7 


j 


bbr + T~ & qqT+T~ 





0.07 ±0.02 


0.06 ±0.02 


0.02 ±0.02 


0.07 ±0 


02 


0.010 


73.7 


k 


qqW & bb^77 


11 


10.2 ±0.7 


5.4 ±0.5 


4.79 ±0.46 


0.34 ±0 


03 


0.219 


73.4 


1 


qqVF & bbvV 





1.46 ± 0.15 


1.23 ±0.13 


0.23 ±0.07 


0.28 ±0 


03 


0.061 


73.4 


m 


qqqq & qqbb 


26 


32.4 ±2.5 


9.1 ±0.6 


23.3 ±2.4 


0.40 ±0 


03 


0.300 


77.7 


n 


qqbb & qqqq 


6 


4.33 ±0.41 


2.73 ±0.18 


1.60 ±0.37 


0.19 ±0 


02 


0.189 


77.7 





qqbb & qqqq 


2 


4.26 ±0.32 


3.16 ±0.19 


1.10 ±0.26 


0.22 ±0 


02 


0.189 


77.7 



Table 7: The 200 GeV data. The description of the entries is given in the caption of Table 2. 
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y/s = 202 GeV 





Selection 


^obs 






^back 


^chan 




Bzz 


Lint 


a 


f+frf+fr 





0.27 ±0.05 


0.22 ±0.01 


0.04 ±0.05 


0.56 ±0.02 


0.010 


38.3 


b 


c^cTvv 





0.37 ±0.04 


0.23 ±0.03 


0.13 ±0.03 


0.46 ±0 


06 


0.013 


37.1 


c 







0.44 ±0.05 


0.23 ±0.03 


0.22 ±0.04 


0.44 ±0 


06 


0.013 


37.1 


d 


qqe + e~ &bbe+e~ 


2 


0.99 ±0.06 


0.85 ±0.04 


0.14 ±0.05 


0.62 ±0 


03 


0.037 


36.6 


e 


qqe + e~ &bbe + e~ 


1 


0.19 ±0.02 


0.17 ±0.02 


0.02 ±0.01 


0.46 ±0 


05 


0.010 


36.6 


f 


qq/i + /i~ & bb/z+/z _ 





1.10 ±0.04 


1.01 ±0.03 


0.09 ±0.03 


0.73 ±0 


02 


0.037 


36.6 


g 


qq/i + /i _ & bb\x + \x~ 





0.25 ± 0.02 


0.24 ±0.01 


0.01 ±0.01 


0.64 ±0 


05 


0.010 


36.6 


h 


qqr + T~ & bbr+r - 


1 


0.47 ±0.04 


0.41 ±0.04 


0.06 ±0.02 


0.30 ±0 


03 


0.037 


36.6 


i 


qqr+r - & bbr+r - 





0.12 ±0.02 


0.11 ±0.02 


0.01 ±0.01 


0.28 ±0 


04 


0.010 


36.6 


j 


bbr + T~ & qqr+T~ 





0.04 ±0.01 


0.03 ±0.01 


0.01 ±0.01 


0.07 ±0 


02 


0.010 


36.6 


k 


qqVF & bbvV 


3 


4.66 ± 0.32 


2.53 ±0.23 


2.13 ±0.23 


0.31 ±0 


03 


0.219 


36.3 


1 


qqvV & bbvV 


1 


0.71 ±0.07 


0.64 ±0.07 


0.08 ±0.03 


0.28 ±0 


03 


0.061 


36.3 


m 


qqqq & qqbb 


22 


21.5 ± 1.7 


5.4 ±0.4 


16.1 ±1.7 


0.48 ±0 


03 


0.300 


36.6 


n 


qqbb & qqqq 


2 


1.54 ±0.14 


1.03 ±0.08 


0.51 ±0.12 


0.15 ±0 


01 


0.189 


36.6 





qqbb & qqqq 





2.21 ±0.17 


1.63 ±0.10 


0.59 ±0.14 


0.23 ±0 


02 


0.189 


36.6 



Table 8: The 202 GeV data. The description of the entries is given in the caption of Table 2. 



y/s = 205 GeV 





Selection 


T^obs 




nzz 


^back 


^chan 




Bzz 




a 


£+£-£+£- 


1 


0.73 ±0.12 


0.49 ±0.01 


0.24 ±0.12 


0.55 ±0 


02 


0.010 


81.7 


b 


c^cTvv 





0.68 ±0.09 


0.47 ±0.06 


0.21 ±0.06 


0.39 ±0 


05 


0.013 


84.4 


c 


H + H~ vT> 





0.81 ±0.11 


0.48 ±0.07 


0.33 ±0.08 


0.40 ±0 


06 


0.013 


84.4 


d 


qqe + e~ &bbe+e~ 


1 


2.19 ±0.12 


1.83 ±0.09 


0.36 ±0.08 


0.59 ±0 


03 


0.037 


80.5 


e 


qqe + e~ &bbe + e~ 


1 


0.51 ±0.05 


0.44 ±0.04 


0.06 ±0.03 


0.51 ±0 


05 


0.010 


80.5 


f 


qq/i + yU~ & bb,u+/!~ 


4 


2.56 ±0.09 


2.43 ±0.08 


0.13 ±0.05 


0.78 ±0 


03 


0.037 


80.5 


g 


qq/i + /i~ & bb/i + /i~ 


2 


0.51 ±0.04 


0.47 ±0.03 


0.04 ±0.03 


0.55 ±0 


05 


0.010 


80.5 


h 


qqT + T~ & bbr+T~ 





1.06 ±0.10 


0.92 ±0.08 


0.13 ±0.05 


0.30 ±0 


03 


0.037 


80.5 


i 


qqr + T~ & bbr + T~ 





0.31 ±0.04 


0.25 ±0.04 


0.06 ±0.02 


0.29 ±0 


05 


0.010 


80.5 


j 


bbr + T~ & qqr+T~ 


1 


0.06 ±0.02 


0.05 ±0.02 


0.01 ±0.01 


0.06 ±0 


02 


0.010 


80.5 


k 


qqW & bb^77 


7 


10.9 ±0.7 


5.8 ±0.5 


5.10 ±0.50 


0.31 ±0 


03 


0.219 


79.9 


1 


qqVF & bbvV 





1.68 ±0.18 


1.48 ±0.15 


0.20 ±0.09 


0.29 ±0 


03 


0.061 


79.9 


m 


qqqq & qqbb 


36 


25.1 ± 1.9 


8.0 ±0.6 


17.0 ± 1.8 


0.32 ±0 


02 


0.300 


80.2 


n 


qqbb & qqqq 


3 


4.46 ±0.42 


2.82 ±0.18 


1.64 ±0.38 


0.18±0 


01 


0.189 


80.2 





qqbb & qqqq 


4 


3.18 ±0.24 


2.51 ±0.17 


0.67 ±0.18 


0.16 ±0 


01 


0.189 


80.2 



Table 9: The 205 GeV data. The description of the entries is given in the caption of Table 2. 
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y/s = 207 GeV 





Selection 


^obs 


n-SM 


nzz 


^back 


^chan 


£?zz 




















(pb- 1 ) 


a 


£+£-£+£- 


2 


1.23 ±0.20 


0.83 ±0.03 


0.40 ±0.20 


0.55 ±0.02 


0.010 


136.9 


b 


c + e~ vv 


1 


1.16 ±0.14 


0.81 ±0.11 


0.36 ±0.10 


0.41 ±0.05 


0.013 


137.6 


c 







1.41 ±0.19 


0.79 ±0.11 


0.62 ±0.15 


0.40 ±0.06 


0.013 


137.6 


d 


qqe + e~ &bbe+e~ 


6 


3.72 ±0.21 


3.12 ±0.16 


0.59 ±0.13 


0.59 ±0.03 


0.037 


134.0 


e 


qqe + e~ & bbe + e~ 





0.82 ±0.08 


0.72 ±0.07 


0.10 ±0.05 


0.49 ±0.05 


0.010 


134.0 


f 


qq/i + /i~ & bb/i+/i~ 


1 


4.37±0.15 


4.16 ±0.13 


0.21 ±0.08 


0.79 ±0.03 


0.037 


134.0 


g 


qq/i + /i~ & bb/i + /i~ 





0.82 ±0.07 


0.75 ±0.05 


0.07 ±0.04 


0.51 ±0.05 


0.010 


134.0 


h 


qqr + r~ & bbr+T~ 


2 


1.83 ±0.16 


1.61 ±0.14 


0.22 ±0.08 


0.30 ±0.03 


0.037 


134.0 


i 


qqr + r _ & bb-r + T~ 





0.49 ±0.08 


0.40 ±0.06 


0.09 ±0.04 


0.27 ±0.05 


0.010 


134.0 


j 


bbr + r~ & qqr+T - 





0.09 ±0.03 


0.08 ±0.03 


0.01 ±0.01 


0.06 ±0.02 


0.010 


134.0 


k 


qqyV & hhvv 


16 


18.9 ± 1.2 


10.1 ±0.9 


8.80 ±0.86 


0.32 ±0.03 


0.219 


133.6 


1 


qqVz7 & bbvv 


7 


2.88 ±0.30 


2.53 ±0.26 


0.35 ±0.15 


0.29 ±0.03 


0.061 


133.6 


m 


qqqq & qqbb 


36 


41.8 ±3.1 


13.6 ±0.9 


28.3 ±3.0 


0.32 ±0.02 


0.300 


133.3 


n 


qqbb & qqqq 


7 


7.48 ±0.69 


4.75 ±0.30 


2.73 ±0.63 


0.18 ±0.01 


0.189 


133.3 





qqbb & qqqq 


6 


5.34 ±0.41 


4.23 ±0.28 


1.11 ±0.29 


0.16 ±0.01 


0.189 


133.3 



Table 10: The 207 GeV data. The description of the entries is given in the caption of Table 2. 
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